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We study the effect of spin injection into s— and d— wave superconductors, with an emphasis on the 
interplay between boundary and bulk spin transport properties. The quantities of interest include 
the amount of non-equilibrium magnetization (m), as well as the induced spin-dependent current 
(7 S ) and boundary voltage (V s ). In general, the Andreev reflection makes each of the three quantities 
depend on a different combination of the boundary and bulk contributions. The situation simplifies 
either for half-metallic ferromagnets or in the strong barrier limit, where both V s and m depend 
solely on the bulk spin transport/relaxation properties. The implications of our results for the 
on-going spin injection experiments in high T c cuprates are discussed. 

PACS numbers: 71.10.Hf, 73.40. -c, 71.27. +a, 72.I5.Gd 



A number of spin-injection experiments have recently 
been carried out for the high T c cuprates Jj]-^]. These 
experiments are of interest for a variety of reasons. In 
particular, they in principle allow us to extract the bulk 
spin transport properties of the high T c cuprates. In the 
normal state, such spin transport properties have been 
proposed as a probe of spin-charge separation |ffln]. In 
the superconducting state, they can provide important 
clues to the nature of the quasiparticlcs. In addition, 
the anisotropy of the spin transport properties should 
shed new light on the nature of the c— axis transport. 
Finally, spin-injection into superconductors also provides 
a setting to study the phenomenon of spatial separation 
of charge and spin currents K|. 

For the purpose of extracting the bulk spin transport 
properties, it is essential that their contributions to the 
measured physical quantities are separated from those of 
the boundary transport processes. The boundary of in- 
terest here is formed between the high T c cuprates and 
ferromagnetic metals. One transport process through 
such a boundary is the Andreev reflection jlc| ], which 
carries pair-current into the superconductor. The novel 
features of the Andreev reflection involving a ferromagnet 
have recently been addressed |ll|-[l^1 . Another boundary 
process is the single-particle transport. While the charge 
transport involves both processes, the spin transport pro- 
ceeds through the single-particle process only. The inter- 
play between the two processes is therefore expected to 
play an important role in the spin injection experiments. 
The consequences of such an interplay are explored in 
this paper. The s— wave case is simpler, which we will 
address first. (This part of our analysis is also relevant 
to the spin injection experiments in the non-cuprate su- 
perconductors |l4|,[l5|].) We will then extend our analysis 
to the d— wave case. 

Several physical quantities are of interest in spin injec- 
tion experiments. One is the amount of non-equilibrium 
magnetization (m) injected into the superconductor. 
The others include the spin-dependent current (7 S ) and 
boundary voltage (Vs) induced by the injected magneti- 



zation. We will show below that, in general these quan- 
tities reflect very different combinations of the bulk spin 
transport and interface transport contributions. 

The relevant experimental setups are illustrated in Fig. 
|l|. Fig. |l|a) is representative of those used in the on- 
going spin-injection experiments in the high T c cuprates 
An injection current (I, per unit area) from a 
ferromagnetic metal (FM1) to the superconductor (S) is 
applied, and the critical current of the superconductor is 
then measured. The suppression of the critical current, 
AJ C ee (J C (I = 0) - J c (f))/J c (I = 0), is expected to be 
a measure of the amount of injected magnetization (m). 
Fig. ||d) illustrates the spin-injection-detection setup of 
Johnson and Silsbee [j^l). Here a superconductor (S) 
is in contact with two itinerant ferromagnets, FM1 and 
FM2. The magnetization of FM1 is either parallel (er =f) 
or antiparallel (a =1) to that of FM2. For a given a, I a is 
the induced current across the S-FM2 interface in a closed 
circuit. Likewise, V a is the induced boundary voltage 
(V<j) in an open circuit. The spin-dependent current and 
boundary voltage are defined as I s = if — 7| and V s — 
V[ — V± , respectively. 

To highlight the interplay between the boundary and 
bulk transport processes, we will make a number of sim- 
plifying assumptions. The superconductor (S) is assumed 
to be a BCS superconductor, with either an s— wave 
or a d x 2_ y 2— wave order parameter. The ferromagnetic 
metal will be simply modeled by an exchange energy h 
0,0. The Hamiltonians for FM1 and FM2 are, Hi = 

J2k,r( £ k ~ hT ) c k,T C k,T and H.2 = J2k.r( e k ~ hTa ) c l,r C k,r, 

respectively. The unshifted energy dispersions for both 
ferromagnets are assumed to be = h 2 k 2 /2m, as is the 
normal state energy dispersion for the superconductor. 
In addition, the Fermi wavevectors of the superconduc- 
tor and the ferromagnets in the absence of polarization 
are assumed to be equal to kf. For FM1, this implies 
kf T — kft T and a Fermi velocity Vf T — Vft T where 



t T — J 1 + rh/tf 



(1) 
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I, 
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FIG. 1. a) A spin injection setup involving a fer- 
romagnetic metal (FM1) and a superconductor (S). b) A 
spin-injection-detection setup involving a superconductor (S) 
and two ferromagnetic metals (FM1 and FM2). 



For FM2, kf T = kft aT and Vf T — Vft aT . 

We introduce m(x) to denote the steady state spin 
magnetization density in the superconductor. Since m{x) 
is entirely carried by the quasiparticles, it satisfies 



D s d 2 m/dx 2 = -m(a;)/Ti 



(2) 



where T\ is the longitudinal spin relaxation time and D s 
is the spin diffusion constant. 

We now specify the boundary conditions. For the setup 
of Fig. la), x — and x — d correspond to the FM1-S 
interace and the open end of S, respectively. For Fig. lb), 
they instead describe the FM1-S and S-FM2 interfaces, 
respectively. At x = 0, 



D s dm/dx\ x = = I. 



spin 



(3) 



For the setup of Fig. la), and to the leading order in Fig. 
lb) as well, the boundary condition at x = d is 



- D s dm/dx\ x=d = 



(4) 



The crucial question then is, what is the spin-current 
Ispin for a fixed injection electrical-current I. At the 
FM1-S boundary the electrical current is the sum of a 
pair current (I pa ir), carried through Andreev reflection, 
and a single-particle current (I S p), 

I Ipair ~t~ 1 sp (5) 

The spin current, on the other hand, has no contribution 
from the Andreev process given that the Cooper pairs 
are spin singlets. We then expect, 



where rj is the single-particle current polarization, which 
itself also needs to be determined. In the following, we 
first calculate the three currents, I pa ir, hp, and I S pin, for 
a given voltage, V, across the FM1-S barrier. The latter 
can then be expressed in terms of /, through Eq. (j^), 
leading to an expression for I sp i n as a function of /. 

For a fixed V, the values of the three currents depend 
on the nature of the interface. Here, we model this in- 
terface by a delta-function potential, TLinterface — VS(x). 
The Bogoliubov-deGennes equation ftiM is, 



TLq — hr 
A* 



-(Ho + hr) 



E 



(7) 



where TLq = p 2 /2m — /i + TLinterface- Here, h and A are 
non-zero only for x < and x > 0, respectively. In- 
side FM1, the solution is the sum of three plane waves, 
describing the incident electron of spin r and wavevec- 
tor ki T , the Andreev-reflected hole of spin f, wavevector 
k aT and amplitude a T , and the normally-reflected elec- 
tron of spin r, wavevector kb T and amplitude b T , respec- 
tively. Inside the superconductor, the solution is the sum 
of two plane waves, one without branch-crossing and of 
wavevector k CT and amplitude c T , the other with branch- 
crossing and of wavevector kd T and amplitude d T . The 
wavevectors are determined by the requirement that, the 
components parallel to the interface are equal. The am- 
plitudes are then calculated by solving Eq. (Q), together 
with the boundary conditions |fl8f appropriate to this 
equation. These amplitudes, in turn, allow us to deter- 
mine the various currents. 

In the linear response regime, the pair, single- 
particle, and spin currents can be written as I\/V = 

fdE(-df/dE)Z^S(C k - E)vti\{E,k). Here, ^ = 
efe — e/ — hr, f(E) is the Fermi-Dirac distribution func- 
tion, J- is & filtering factor which specifies the distribution 
of incident angles. In addition, 



i T — Op 1 A 
l pair ^ ^-t 



il= e 2 (C T + D T ) 



(8) 



where C T = 
\d T \ 2 {u{k dT ) 2 - 



\c T \ 2 (u(k CT ) 2 - v(k CT ) 2 )v]/v x kT , D T = 
v{k dT ) 2 )v x f /v x kT , and A T = \a T \ 2 v x kf / v x kT . 
(Here kr labels the wavevectors satisfying = E.) 

The potential barrier can be represented in terms of 
a dimensionless quantity, z = m V/h 2 k f . In the follow- 
ing, we will focus on the two extreme limits, z —> and 
z — > oo. In the latter case we expand in terms of l/z. 
In addition, we will consider only temperatures low com- 
pared to the superconducting gap A. These cases are suf- 
ficient to illustrate the main points of this paper. More 
general cases will be discussed elsewhere. 
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s—wave, vanishing barrier: We address first the case 
of an isotropic s—wave superconductor, in the limit of a 
vanishing barrier, z = 0. For k B T <C A, we found, 



Ipair/V— goPsfi 

I sp /V= 9oS s fiY sfi (T) 
I sp in/V= [pb B /e)goM s>0 Y sfi (T) 



(9) 



Here go = Novje 2 where No is the density of states at the 
unshifted Fermi energy, P s ,o = 16a pa i r titi/(t^ti + l) 2 , 
S s , = 2^a sp [t } {t 2 + 1) + t L (t 2 + l)]/(t T *4. + !) 2 > and 
M s , - 2^a spm m -t i )(l-t T t i )/(t T t i + l) 2 . The 
(polarization-dependent) factors a pa ir, and a sp i n are 
of order unity; they describe the averaging over the an- 
gle of incident electrons, and are all equal to unity when 
forward incidence dominates. 

Both I sp and I sp i n have an exponential temperature 
dependence, as specified by the function Y St o(T) = 
v / k B T/Ae~ A/kBT . This reflects the simple physics that 
the transport of both the single-particle current and spin 
current involve the quasiparticles of the superconductor 
and hence only energies above the superconducting gap. 

We are now in a position to determine the non- 
equilibrium magnetization (m) and the induced spin- 
dependent boundary voltage (V s ) and current (J s ) as 
a function of /. Consider first the general case, when 
the ferromagnet has a less than 100% polarization. For 
k B T <C A, Ipair then dominates the total current in Eq. 
(||). As a result, 

hpin/I = {(i B /e)(M S! o/P s , )Y sfi (T) (10) 

Combining Eq. @ with Eqs. (|, |, |), we have, 

m(x = d)/I = fa B /e)(M afi /P s , )Y s ,o(T) . T ) (11 

o s smn(a/o s ) 



where S s = ^JD S T\ is the spin-diffusion length of the 
superconductor. 

This magnetization accumulation leads to a drop 
across the interface of an effective magnetic field [ pi] , 
AH = m/x where x is the uniform spin susceptibility of 
the superconductor. In the spin-injection-detection setup 
(Fig. lb), for a closed circuit, this effective field drop will 
induce a current across the S-FM2 interface. Following a 
procedure similar to the one leading to Eq. (ph, we found 



rclosed 



= vgo(VB/e)M sfi Y S}0 (T)(m/x) 



(12) 



The induced spin-dependent current, I s = I^ losed 
jciosed^ j g t ken equal to 



I s /I - 2 ffo (MB/e) 2 (M s yP s ,o)[F s ,o(T)] 



current that would have been induced by the effective 
magnetic field drop. The total induced current is, 

I° pen = ago^B/e)M s ,oY sfi (T)(m/ X ) + g P s , Q V a (14) 



Compared to Eq. (|12|), the additional term is the current 
induced by the voltage drop; here only the pair current 
has been kept as it dominates over the corresponding 
single-particle current. Setting I°P en = determines V a , 
leading to a spin-dependent boundary voltage, V s = V[ — 
Vi, as follows, 



V s /I = 2(p B /e) 2 (M sfi /P s ,o) 2 [Y s ,o(T)f 



Ti 



(15) 



X§ s sinh(d/(5, 

Eqs. (0, H ||) reveal one key point: In general each 
of the three quantities, m, V^, and I s , reflects a different 
combination of the boundary and bulk contributions. 

Consider next the special case of a half-metallic ferro- 
magnet. Ipair vanishes in this case, as can be seen from 
the expression for P Si o in Eq. (^); the Andreev reflection 
is completely suppressed in this case |lM. The total cur- 
rent in Eq. (|5|) is then entirely given by I sp . This, com- 
bined with the fact that M S: o = >5 Sl o for a half-metallic 
ferromagnet, leads to a simple relationship between the 
spin and charge current, 

hpin/I ' = fi B /e (16) 

In addition, the second term in Eq. (0) is now replaced 
by the corresponding single-particle current. The results 
then become, 



m(x = d)/I= (n B /e) 



5 S sinh(<i/(5 s 



I s /I= 2 90 {fi B /e) 2 M Si0 Y sfi (T) 



Ti 



X$ s sinh(d/(5 s ) 



V s /I= 2(/i fl /e) 1 



xS s Biah(d/6 a ) 



(17) 



xS s sinh(d/^ s ) 



(13) 



Similarly, in the open circuit case, a boundary voltage 
will develop across the S-FM2 interface to balance the 



Eq. ( |l7|) reveals another main conclusion of this paper: 
The boundary-transport-contributions are absent in the 
expressions for the non-equilibrium magnetization, m, 
and the induced spin-dependent boundary voltage, V s , 
when the ferromagnet is half- metallic. For m, this is the 
result of the simple relationship between the spin current 
and total current (Eq. (|H|)), which arises whenever the 
Andreev reflection is absent. (Note that the total current 
I is fixed.) For V s , the reasoning leading to this conclu- 
sion is slightly more subtle. When the Andreev reflection 
is absent, the boundary-transport-contributions give rise 
to exactly the same prefactors in the two currents in Eq. 
(|l4|). Since the boundary voltage V a is determined by 
balancing these two terms, the boundary-transport fac- 
tors cancel out exactly in the expression for V s . 

The above cancellation argument does not apply to 
the induced current I a . Indeed, / s does contain the 
boundary-transport factors. 
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s—wave, large barrier: Consider now the limit of 
a large barrier. The results for the pair, single-particle, 
and spin currents parallel Eq. (^), except that the factors 
P sfi , S S: o, M S: o, and Y Si0 are replaced by 

P s ,oc = (1/z )5oQ!poir*T*J. 

V2(l/z 2 ) 50 a sp (i T +t l ) 



(18) 




/ n/2(l/z )g a sp m 



respectively. The resulting expressions for m, I s and 
are also given by Eqs. ( pT| , |i~3| , [j~5| ), with an appropriate 
substitution by the factors given in Eq. (p^|). 

The pair current is of order 1/z 4 . Both the single- 
particle current and spin current, on the other hand, are 
of order 1/z 2 . The latter is to be expected, since in the 
large barrier limit the single-particle transport can be 
described in terms of a tunneling picture. For the same 
reason, the temperature dependence of I sp and I sp m in 
this case should reflect simply the thermal smearing of 
the single-particle density of states in the superconduc- 
tor. Indeed, Y St00 (T) corresponds to the low temperature 
limit of the Yosida function. 

Whenever the Andreev reflection is non-negligible, 
each of the three quantities, m, V s and J s , again depends 
on a different combination of the boundary transport and 
bulk transport properties. 

When the potential barrier is so strong that the 
temperature-independent 1/z 4 contribution is negligible 
compared to the temperature-dependent 1/z 2 terms, the 
Andreev process is absent. Here again, while the bound- 
ary transport terms still appear in J s , they are cancelled 
out in to and V s . 

d—wave case: We now turn to the case of a d x 2_ y 2 
superconductor. Here we will focus on the case when 
both the FM1-S and S-FM2 interfaces involve the {110} 
surface of the cuprates. Consider first the limit of a van- 
ishing barrier. The result for the currents again parallels 
Eq. (^), with the prefactors replaced, respectively, by 

Pd,0= ISgoCtpairt^ti/fa +t i ) 2 

S d , = (ln2)g a sp [C^ + (T— 4-)]/(*T*4. + l f 
M d<0 = (ln2)g a spm [C^ - (W)]/(Ml + I) 2 
Yd,o{T)= k B T/A (19) 



where C T = (t l + l) 2 + 4(3 - tt)^ - (19/3 - 2tt)(^ 
For the large barrier limit, we found Pd.oo = 
Y dt00 (T) = Y dt0 (T), and 



M d , 



(4-7r)(ln2/2)(l/z 2 )a sp (i T +t i ) 
(4-7r)(ln2/2)(l/z 2 )a spm (i T - t J 



-I) 2 . 
Pd,o, 



(20) 



To summarize, we have studied the effect of spin in- 
jection into s— and d— wave superconductors. Through 
the explicit results in the small and larger barrier lim- 
its, we conclude that Andreev-reflection makes the differ- 
ent physical properties measured in spin injection experi- 
ments reflect very different combinations of the boundary 
and bulk spin transport properties. For the purpose of 
isolating bulk spin transport properties, it is desirable 
to make the Andreev contribution as small as possible. 
In this case, the non-equilibrium magnetization and the 
spin-dependent boundary voltage depends solely on the 
bulk spin transport properties. This occurs in two limits. 
One is the limit of a strong barrier in tunneling geome- 
tries such that the Andreev bound states are absent. The 
other is for the half-metallic ferromagnet. 

This work has been supported in part by the NSF 
Grant No. DMR-9712626, a Robert A. Welch grant, and 
an A. P. Sloan Fellowship. While this paper was under 
preparation we learned of a preprint p0[| of Kashiwaya 
et al. who calculated the non-linear spin current for the 
case of a d— wave superconductor. 



Note that, while I sp and I sp i n 



are still of order 1/z 2 , I pair 
is of order l/z°. The latter reflects the formation of the 
Andreev-bound states [fl9|-[22|. 



[9 

[10 
[11 

[12 
[13 
[14 
[15 

[ie; 
[17; 
[is; 
[19; 

[20 

[21 
[22 
[23 



V. A. Vas'ko et al., Phys. Rev. Lett. 78, 1134 (1997). 

Z. W. Dong et al., Appl. Phys. Lett. 71, 1718 (1997). 

V. A. Vas'ko et al, Appl. Phys. Lett. 73, 844 (1998). 

N.-C. Yen et al., preprint (1998). 

N. Hass et al., Physica C235-240, 1905 (1994). 

Q. Si, Phys Rev. Lett. 78, 1767 (1997). 

Q. Si, Phys. Rev. Lett. 81, 3191 (1998). 

H. L. Zhao and S. Hershfield, Phys. Rev. B52, 3632 

(1995). 

S. A. Kivelson and D. Rokhsar, Phys. Rev. B41, 11693 
(1990). 

A. F. Andreev, Sov. Phys. JETP 19 1228 (1964). 

M. J. M. de Jong and C. W. J. Beenakker, Phys. Rev. 

Lett. 74 1657 (1995). 

J.-X. Zhu, B. Friedman, and C. S. Ting, pre print(1998). 



I. Zutic and O. T. Vails, |cond-mat/9808285 . 
R. J. Soulen et al., Science 282, 85 (1998). 
S. K. Upadhyay et al., Phys. Rev. Lett. 81, 3247 (1998). 
M. Johnson, Appl. Phys. Lett. 65, 1460 (1994) and ref- 
erences therein. 

The spin waves have a super-ohmic spectral function and 

are hence unimportant for our analysis. 

G. E. Blonder, M. Tinkham, and T. M. Klapwijk, Phys 

Rev. B 25, 4515 (1982) and references therein. 

Chia-Ren Hu, Phys. Rev. Lett. 72, 1526 (1994). 

Y. Tanaka and S. Kashiwaya, Phys Rev. Lett., 74, 3451 

(1995); S. Kashiwaya et al., Phys Rev. B 53, 2667 (1996). 

J. H. Xu et al., Phys. Rev. B 53, 3604 (1996). 

M. Fogelstrom et al., Phys. Rev. Lett. 79, 281 (1997). 

S. Kashiwaya, Y. T anaka, N. Yoshida, and M. R. Beasley, 

cond-mat/9812160. 



4 



